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ABSTRACT 


An  experiment  has  been  designed  to  test  the  basic  assumptions  of 
a  new  model^  which  describes  debris  motion  towards  the  northern  conjugate 
region  from  high  altitude  nuclear  explosions.  The  experimental  facility 
.consists  of  a  coaxial  plasma  source  which  accelerates  ionized  aluminum 
into  an  environment  where  the  background  pressure  and  magnetic  field 
can  be  controlled.  The  initial  plasma  velocity  was  10  cm/fisec,  the  back¬ 
ground  magnetic  field  was  varied  from  zero  to  220  gauss  and  the  pressure 
between  2  /i  and  50jiin  air.  These  conditions  are  such  that  (1)  ordinary 
momentum  transfer  collisions  are  unimportant  (2)  the  magnetic  pressure 
is  negligible  as  compared  to  dynamic  plasma  pressure  and  (3)  the  only 
momentum  exchange  mechanism  is  charge  exchange  collisions  in  the 
presence  of  a  magnetic  field.  Ultraviolet  radiation  from  the  plasma  source 
preionizes  the  background  gas  producing  an  Alfven  speed  based  on  the  ion 
density  which  is  less  than  the  plasma  velocity.  The  ambient  magnetic  field 
is  swept  up  by  aluminum  plasma  and  provides  the  mechanism  for  charge 
exchange  momentum  exchange  process.  We  find  that  the  aluminum  plasma 
is  decelerated  with  the  addition  of  a  magnetic  field  nearly  parallel  to  the 
plasma  motion.  The  deceleration  rate  is  proportional  to  the  background 
pressure.  Finally,  we  note  that  the  charge  exchange  pickup  is  effective 
even  at  relatively  low  magnetic  fields  where  the  air -ion  gyro  radius  based 
upon  initial  plasma  velocity  is  several  times  the  plasma  scale  length. 
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I.  INTRODUCTION 


In  the  study  of  high  altitude  nuclear  explosions,  the  motion  of  ionized 
debris  along  the  earth’s  magnetic  field  is  not  well  understood.  It  is  vital 
to  the  problem  of  radar  blackout  that  we  be  able  to  predict  the  debris  deposi¬ 
tion  resulting  from  such  explosions. 

Several  theories  and  calculations  have  been  constructed  to  explain 
this  motion  and  it  is  the  purpose  of  this  experimental  program  to  simulate 
some  of  the  effects  associated  with  high  altitude  nuclear  explosions  in  order 
that  the  basic  assumptions  of  these  theories  may  be  better  understood.  To 
this  end,  we  produce  an  unsteady,  high  velocity  flow  of  ionized  material 
along  a  static  magnetic  field  in  an  environment  where  ordinary  momentum 
transport  collisions  are  unimportant.  In  this  situation,  transfer  of  momentum 
can  only  occur  as  a  result  of  charge  exchange  and  ionizing  collisions  in  the 
presence  of  a  magnetic  field. 

The  focal  point  of  this  experiment  is  the  use  of  electromagnetic 
forces  in  the  form  of  a  magnetic  accelerator  to  ionize  and  accelerate  a 
known  amount  of  metallic  ion  plasma.  With  this  type  of  plasma  accelerator 
it  is  possible  to  impart  many  thousands  of  joules  to  milligram  masses, 
thereby  attaining  speeds  relevant  to  the  study  of  high  altitude  nuclear  ex¬ 
plosions.  Thi3  capability  is  not  presently  available  with  other  energy 
sources  such  as  converging  laser  beams. 

To  perform  a  relevant  laboratory  experiment,  the  following 
criteria  are  necessary: 

„  1.  The  velocity  of  the  ionized  material  must  be  of  the  order  of 

10'  cm/sec  or  higher.  At  this  velocity,  the  charge  exchange  and  direct 
ionization  cross  sections  are  much  greater  than  those  of  ordinary  gas 
dynamic  collisions, 

2.  The  background  gas  pressure  must  lie  in  the  range  1  p  <  p  <  45  p. 
The  upper  bound  insures  that  the  mean  free  path  for  gas  dynamic  collisions 
is  much  larger  than  the  flow  scale  (30  cm).  The  lower  bound  is  added  to 
provide  sufficient  interaction  with  the  background  gas  within  the  scale  of 
the  experiment  (-  3  meters). 

3.  The  ion  gyro  radius  must  be  smaller  than  or  equal  to  the  flow 
scale  so  that  the  magnetic  field  can  play  an  effective  role  in  the  momentum 
transfer  process. 

4.  The  magnetic  energy  of  the  static  field  must  be  much  less  than 
the  kinetic  energy  of  the  plasma  so  that  the  magnetic  field  does  not  directly 


influence  the  plasma  motion.  The  magnetic  field  is  solely  an  intermediary 
agent  in  the  momentum  transfer  process  as  is  the  case  in  high  altitude 
nuclear  debris  motion. 


5.  The  energy  associated  with  the  convected  B  field  imbedded  in 
IJjmetal  plasma  must  also  be  smaller  than  the  kinetic  energy  of  the  debris. 
‘'We  wish  to  study  the  motion  of  ionized  debris  without  the  added  complica¬ 
tion  of  a  "magnetic  piston"  behind  it.  The  electric  currents  in  the  plasma 
necessary  to  sustain  the  jump  in  magnetic  field  would  unduly  complicate 
*  the  experiment. 


'  6.  It  is  desirable  to  have  a  preionization  of  the  ambient  air  available 

:  before  the  metal  ion  plasma  interacts  with  it.  This  will  establish  a  relatively 
T  low  Alfven  speed  ahead  of  the  debris,  based  on  the  ion  density  and  3  field 
L’;  ahead  of  the  expansion.  The  velocity  of  the  metal  ion  plasma  is  greater  than 
the  Alfven  speed  ahead,  similar  to  the  conditions  for  actual  high  altitude 
debris  motion.  In  this  experiment,  the  violent  discharge  in  the  accelerator 
creates  sufficient  ultra-violet  radiation  to  ionize  five  percent  of  the  back¬ 
ground  gas  at  a  distance  of  1  meter  from  the  plasma  source. 


The  most  significant  question  to  be  answered  by  this  experiment  is: 
Can  the  magnetic  field  which  is  aligned  parallel  to  the  flow  affect  the  motion 
of  the  ionized  debris  and  if  so,  how,  in  detail,  is  this  accomplished?  The 
second  half  of  this  question  involves  a  knowledge  of  the  magnetic  field  con¬ 
figuration  which  accompanies  the  debris  motion.  Once  this  magnetic  field 
interaction  is  established,  it  then  remains  to  determine  the  functional 
behavior  upon  the  parameters  of  background  gas  pressure  and  ambient 
magnetic  field. 


II.  DESIGN  OF  EXPERIMENT 


The  power  requirements  for  the  plasma  accelerator  are  essentially 
fixed  by  the  plasma  velocity  and  the  requirement  that  the  ion  gyro  radius 
be  smaller  than  or  comparable  to  the  flow  scale  in  the  following  manner. 
The  power  unit  area  of  the  expanding  debris  bubble  is  given  by  the  dynamic 
pressure  of  the  ambient  gas  multiplied  by  the  expansion  velocity.  Thus, 
the  total  power  expended  by  the  expanding  bubble  is  approximately 

P  *  (pv2)  (v)  (4;rR2)  (1) 

where  p  =  background  density 
v  =  expansion  velocity 
R  =  bubble  radius 

The  background  density  is  fixed  by  restriction  on  the  ion  gyro 
radius,  i.e.,  that 


r . 


(2) 


where  r^  is  given  by 


r . 


U) 


(3) 


and  qj  and  m^  are  the  charge  and  mass  of  an  ion.  Also  we  must  have  the 
dynamic  pressure,  pv^,  larger  than  the  ambient  magnetic  field  pressure 
b^/2Mq-  i*  e-  « 


2 


Mo 


(4) 


Solving  for  p  and  substituting  into  Eq.  (1),  one  obtains  the  power  for  the 
accelerator 


(5) 
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Assuming  an  average  ion  mass  of  20  and  velocity  of  10  m/sec,  it 
can  be  seen  from  Eq.  (4)  that  the  minimum  (i.e.  ,  R/ri  =  1)  power  require¬ 
ment  is  5  x  10'  p  watts. 

apjljlS  The  time  scale  of  the  experiment  and  therefore  the  energy  input,  is 

?'  .xhosen  to  give  a  debris  bubble  of  convenient  size.  To  expand  at  10?  cm/sec 
to  a  diameter  of  30  cm,  3  ji.  sec  arc  required.  If  one  chooses  (R/r^)2  p  ^  30 
_  then  this  corresponds  to  an  input  energy  of  1500  joules. 

The  design  for  the  plasma  source  which  fulfills  these  power  and 
energy  requirements  is  shown  schematically  in  Fig.  1,  The  power  delivered 
tp  the  accelerator  is 

P  =  VI  (6) 

where  V  voltage  across  the  accelerator  and  I  is  the  current  through  it. 

The  voltage  is  given  by 


V  =  v  B 


TOT  r 


(V) 


where  BTqT 


B0  =  V/U  »  r) 

r  =  anode  radius 


The  z  component  of  the  magnetic  field  which  is  externally  applied  by  a 
coil  surrounding  the  accelerator,  serves  two  purposes;  the  axial  lagnetic 
field  (1)  limits  the  heat  transfer  to  the  anode  and  increases  the  efficiency 
of  the  arc  and  (2)  provides  a  smooth  transition  to  the  ambient  magnetic 
'  T  field  which  is  applied  within  the  test  section. 

r  Under  the  normal  operating  conditions  when  Bz  in  the  accelerator  is 

comparable  to  or  less  than  Bg,  a  power  level  of  10?  watts  corresponds  to 
a  voltage  and  current  of  10  kV  and  100  kA  respectively.  These  input  re¬ 
quirements  are  supplied  by  a  fast  (1  sec  rise  time)  capacitor  bank  con¬ 
taining  10,  000  joules  at  50  kV.  It  follows  that  approximately 


10  kV 
50  kV  x 


10,  000  joules  =  2000  joules 


are  delivered  to  the  plasma  accelerator. 


It  is  desirable  to  have  a  reproducible  mass  of  ions  to  be  accelerated. 
The  amount  of  mass  is  fixed  by  the  input  energy  and  exhaust  velocity  re¬ 
quirements,  i.e., 


2 


v 


2  •  2000  joules 

inl°  2/  2 

10  m  /sec 


.  4  milligram 
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Fig.  1  Plasma  Accelerator 
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and  a  simple  method  of  introducing  this  mass  (in  the  form  of  thin  wires) 
into  the  accelerator  is  shown  in  Fig.  1.  A  better  method  might  be  the 
use  of  foils;  however,  the  low  mass  requirement  forces  the  use  of  very 
j^thin  foils  and  there  is  no  simple  method  of  placing  a  foil  in  the  arc.  The 
jywires,  on  the  other  hand,  are  simply  tied  around  the  cathode  and  rest 
Jupon  the  anode. 

There  are  two  factors  which  aid  in  rapidly  distributing  the  ionized 
material  symmetrically  around  the  arc  prior  to  acceleration.  Both  mech¬ 
anisms  rely  upon  the  observation  that  since  the  propelling  force  on  the 
ionized  aluminum  is  IB q  and  B q  °c  I,  this  force  is  proportional  to  1^. 
Furthermore,  since  I  sin  iot,  at  early  times, 


Force  <x 


where  x  represents  the  axial  position  of  the  ionized  mass.  Upon  perform¬ 
ing  two  quadratures,  it  follows  that  x  t^  and  a  simple  calculation  shows 
that  the  ionized  mass  has  traveled  only  the  order  of  a  few  cm,  by  the  time 
the  current  pulse  has  passed  its  maximum  and  decreased  substantially  to 
zero.  The  ionized  mass  has  not  yet  left  the  vicinity  of  the  accelerator 
when  the  arc  currents  vanish.  The  second  half  cycle  currents  restrike  on 
the  insulator  and  the  plasma  has  kinetic  energy  with  essentially  no  B  field 
energy. 

The  first  symmetrizing  factor  is  the  axial  bias  field.  The  axial 
bias  field  is  supplied  from  a  slow  (5  millisecond)  capacitor  bank,  and  as 
a  result,  this  field  is  essentially  constant  over  the  entire  time  scale  of 
the  discharge.  Now,  at  the  beginning  of  the  arc  discharge  cycle,  when 
the  current  and  therefore,  Bg,  is  low,  the  ionized  material  is  spun  in  a 
tight  spiral  by  the  axial  bias  field.  By  the  time  the  arc  current  is  suf¬ 
ficiently  strong  to  propel  the  ionized  mass  out  of  the  arc  region,  azimuthal 
symmetry  should  be  achieved. 

The  second  factor  aiding  symmetry  is  the  instabilities  associated 
with  exploding  wires.  With  the  passage  of  105  amps,  the  wire  should 
definitely  be  kink  unstable  in  times  shorter  than  1  /i,  sec.  Instabilities  of 
this  type  should  help  distribute  the  plasma  around  the  arc. 

Image  converter  photograp]  s  taken  at  l/2  /x  sec  after  initiation  of 
the  arc  current  with  exposure  times  of  50  nanoseconds  do  not  show  any 
gross  asymmetries  in  the  arc  discharge  and  we  conclude  that  these  wires 
represent  a  simple,  adequate  method  of  reproducibly  introducing  a  known 
amount  of  material  into  the  arc. 

The  pressure  of  the  background  gas  is  determined  in  such  a  manner 
so  as  to  slow  the  debris  in  a  distance  which  is  convenient  with  the  experi¬ 
mental  facility.  The  debris  will  slow  down  appreciably  if  it  picks  up  a  mass 
of  air  equal  to  its  initial  weight  (.  4  mg).  Assuming  a  bubble  diameter  of 
30  cm  and  complete  mass  pickup,  this  pickup  requirements  corresponds  to 
an  air  pressure  of  ~  10(1. 


* 


Figure  2  shows  a  schematic  diagram  of  the  experimental  facility 
including  some  of  the  diagnostic  devices.  The  time  constant  of  the  slow 
capacitor  bank  is  long  enough  to  permit  the  (axial  bias)  magnetic  field 
to  penetrate  the  copper  anode  and  connect  smoothly  with  the  (ambient) 
magnetic  field  produced  by  the  field  coils.  A  time  delay  generator  switches 
on  the  fas  •  bank  when  the  ambient  field  has  reached  its  maximum, 

The  main  diagnostic  devices  are  a  set  of  rotating  mirror  cameras 
and  various  magnetic  field  probes.  The  cameras  have  a  maximum  writing 
speed  of  1  cm/sec  can  resolve  velocities  equal  to  10  cm / \i  sec  and  the  lo¬ 
cation  of  the  two  cameras  provide  an  almost  continuous  field  of  view  for  the 
first  meter  of  travel.  The  magnetic  field  probes  are  "center  tap"  grounded 
and  are  operated  "push-pull"  to  eliminate  electrostatic  pickup. 

An  important  monitoring  device  is  a  tungsten  U.  V.  detector  which 
is  located  on  the  rear  wall  of  the  tank  opposite  the  accelerator.  This 
detector  gives  the  li;vel  of  U.  V.  radiation  resulting  from  the  violent  break¬ 
down  of  the  aluminum  plasma,  in  the  wavelength  range  200-1200.$..  From 
this  information,  the  initial  ionization  level  of  the  background  gas  can  be 
calculated.  It  is  to  be  noted  that  high  altitude  nuclear  explosions  also  pro¬ 
duce  a  high  level  of  background  ionization  as  a  result  of  x-radiation. 
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III.  EXPERIMENTAL  RESULTS 


1.  Accelerator  Performance 

Voltage  and  current  measurements  made  at  the  plasma  accelerator 
indicate  its  performance  is  in  accordance  with  design  specifications.  The 
peak  voltage  across  the  accelerator  was  measured  to  be  ~  5  kV.  The 
current  signal,  on  the  other  hand,  exhibits  a  quarter  cycle  of  T  =  1.5  (l  sec 
which  indicates  a  peak  current  of 

C  V  7T  (. 

I  =?  — — —  =?  2x10  amps 

where  VQ  =  30  kV  is  the  capacitor  voltage.  Thus,  the  average  (RMS)  power 
is  approximately  l/2  x  10^  watts  which  gives  an  input  energy  of  1,  500  joules 
for  the  first  l/2  cycle. 

Assuming  that  all  this  input  energy  is  converted  into  kinetic  energy 
of  the  aluminum  debris  (the  validity  of  this  assumption  has  been  confirmed 
by  optical  and  field  probe  data  and  will  be  discussed  shortly),  the  exhaust 
velocity  of  the  ionized  aluminum  is  approximately 

v  ,  .  ^F- .  *=  10^  m/sec 

exhaust  y  m  ^  input 

.7 

for  m  =  4  x  10  kg. 

The  E/B  drift  velocity  in  the  nozzle  region  of  the  anode  is  consistent 
with  this  exhaust  velocity.  The  azimuthal  magnetic  field  in  the  annulus  is 
approximately 

Ml  /2 

B„  *=  — =5-  ~  1.6  web/m 

o  7 rK 

where  R  =  .  05  m  is  the  anode  radius  and  the  radial  electric  field  is 
approximately 


E  ?  -  ?  2x  105  volts/m 
r  r 

where  r  ="  2.  5  cm  is  the  annular  spacing  and  V  =  5  kV  is  the  accelerator 
voltage.  These  two  orthogonal  fields  given  an  E/B  drift  velocity  in  the 
axial  direction  of  ~  10^  m/sec. 


A  further  check  of  the  accelerator  performance  was  made  by  firing 
the  plasma  into  a  high  pressure  background  with  no  applied  field  and  observ¬ 
ing  the  resultant  motion.  For  background  pressures  of  the  order  of  1  mm, 
•the  density  is  sufficiently  high  so  that  all  of  the  ambient  gas  is  picked  up 
via  standard  momentum  transfer  collisions  and  a  simple  theory  which  con¬ 
serves  momentum  should  predict  the  debris  motion. 

Figure  3  shows  the  data  which  is  an  x-t  diagram  obtained  from  a 
rotating  mirror  camera  along  with  a  pencil  sketch  of  the  data. 

The  philosophy  of  the  calibration  measurement  is  to  calculate  the 
initial  mass  in  the  accelerator  from  the  high  pressure  mirror  camera 
jdata  and  compare  that  with  the  mass  of  aluminum  wire  strung  across  the 

gap- 


The  conservation  of  momentum  when  applied  to  the  aluminum  mass 
when  thermal  effects  are  neglected  and  one  dimensional  flow  is  assumed 
is  given  by 


f  m  +  m(z)  1  v (z )  =  m  v  (8) 

1  O  J  o  o 

where  m  and  v  are  the  mass  and  velocity  and  the  subscript  zero  refers  to 
their  initial  values.  m{z)  is  the  mass  of  ambient  gas  picked  up  by  the 
aluminum  and  is  proportional  to  the  volume  swept  out.  The  volume  swept 
out  at  high  pressures  is  essentially  a  cone  with  base  radius  approximately 
one  third  times  that  of  its  height.  This  data  was  determined  using  an 
S,  T.  L.  image  converter  camera.  The  picked  up  mass  is 

m(z)  =  pair  it  z2/27  (9) 

where  p  .  is  the  air  density  at  1  mm.  Substituting  Eq.  (9)  into  Eq.  (8) 
and  integrating  the  result  once,  the  initial  mass  is  expressed  as 

4 

mo  108  ^air  (v  t-z)  ^ ^ 

o 

The  initial  velocity  is  estimated  from  the  first  10  cm  of  travel  and  the 
resultant  mQ  is  plotted  in  Fig.  3.  The  aluminum  mass  put  into  the  ac¬ 
celerator  is  4  x  10"^  kg.  The  inferred  mass  from  the  mirror  camera 
data  is  fairly  constant  and  agrees  well  for  the  first  20  cm  of  travel.  After 
that,  the  plasma  has  slowed  sufficiently  so  that  its  kinetic  energy  has  been 
converted  to  thermal  energy.  Beyond  this  point,  a  blast  wave  solution 
should  be  used.  In  the  blast  wave  model  plasma  pressure  exert  a  force  on 
the  combined  aluminum  and  air  mass  and  this  decreases  the  deceleration 
and  appears  as  an  effective  increase  in  ir^-,  in  Eq.  (10). 

With  the  simple  method  of  introducing  a  known  amount  of  aluminum 
and  this  calibration  technique,  this  accelerator  can  be  used  as  reliable 
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tool  for  high  velocity  plasma  research.  The  attainment  of  exhaust  velocities 
of  10  m/sec  as  indicated  by  accelerator  voltage  and  current  measurements 
_and  confirmed  by  optical  and  field  probe  data  (in  the  next  section)  show  that 
^the  accelerator  operates  at  high  efficiency  and  that  the  performance  can  be 
predicted  by  simple  energy  and  momentum  considerations. 

2.  Initial  Conditions 

The  electrical  breakdown  of  the  aluminum  wire  along  with  the  passage 
of  the  order  of  10-  amperes  in  each  wire  results  in  the  emission  of  a  pulse 
of  radiation  which  preionizes  a  fraction  of  the  ambient  gas  in  much  the  same 
manner  as  x-radiation  ionizes  the  atmosphere  in  high  altitude  nuclear  ex¬ 
plosions.  A  tungsten  ultraviolet  detector  was  mounted  at  the  rear  of  the 
tank  facing  the  accelerator  and  the  amount  of  radiation  in  the  range  200- 
1200  A  was  recorded. 

The  power  incident  on  the  U.  V.  detector  is  related  to  its  output 
current  in  the  following  manner: 

1 

PINC  =  e  t!x 

where  =  incident  power 

I  =  output  current  (photo  current  from  tungsten  photo  cathode) 

e  =  electronic  charge 

=  photon  efficiency  =  number  of  electrons  emitted  per 
incident  photon 

=  photon  energy 


Assuming  an  average  photon  energy  (in  the  range  200- 1200  J£)  of  15  eV  and 
a  photon  efficiency  of  about  10%,  1  the  incident  power  to  the  detector  was 
found  to  be  ~  30  watts.  The  detector  area  is  A^  =  10"  *  cm  and  it  is  located 
R  =  ~  3m  from  the  accelerator.  Assuming  that  the  power  is  radiated  in  a 
solid  angle  of  27rR  ,  the  peak  power  radiated  by  the  accelerator  in  the 
U.  V.  range  is 


P  =  (30  watts) 


2?rR 


8 


“=  2x10  watts 


which  is  approximately  10%  of  the  input  power  to  the  accelerator. 


To  calculate  the  fractional  ionization  of  the  ambient  air  in  the  tank 
we  note  that  for  photon  energies  greater  than  the  ionization  potential,  the 
photoionization  cross  section,  a,  is  cr  =  10“17  cm^.  Thus  at  pressures 
of  the  orderjDf  lOgt  (density  =  3  x  1014  cm"3),  the  photoionization  mean 
free  path  X  ~  3  meters  which  is  larger  than  the  scale  of  the  experiment. 


Using  a  linear  approximation,  the  fractional  ionization  per  second  is  then 
given  by 


%  ionization/sec  =  (photon  flux)  a 

Since  the  photon  flux  varies  inversely  as  the  square  of  the  distance  from 
the  source,  the  ionization  fraction  becomes  greater  as  one  approaches  the 
accelerator.  The  photon  flux  at  the  U.  V.  detector  is  ~  1.2  x  10^9 
photons/(cm^sec)  and  the  photon  flux  at  any  point  is 

photon  flux  =  1.  2  x  10+^  (— ^  ^  P-— — 

R  cm  sec 

where  R  is  the  distance  to  accelerator  in  meters.  Thus,  since  the  cur¬ 
rents  and  also  the  U.  V.  emission  lasts  for  approximately  10-b  sec,  the 
ionization  fraction  is  the  order  of  5%  at  1  meter  from  the  source  and  100% 
at  20  cm. 

This  calculation  indicates  the  following:  (1)  close  to  the  accelerator 
(~  20  cm)  the  ambient  gas  is  fully  ionized  and  the  pickup  mechanism  must 
be  similar  to  Longmire's  model,  2  (2)  further  from  the  accelerator  (~  lm) 

the  gas  is  95%  neutral  and  the  pickup  mechanism  must  be  associated  with 
charge  exchange  and  ionizing  collisions  in  the  presence  of  a  magnetic  field. 

In  regard  to  the  second  conclusions,  we  note  that  the  cross  section 
for  charge  exchange  and  direct  ionization  collisions  are  of  the  order  10"  15 
cm^  for  ion  energies  of  the  order  1.  5  kev.  The  mean  free  path  for  these 
processes  at  10/i  is  ~  3  cm  which  is  much  smaller  than  the  flow  scale. 

The  mean  free  path  for  ordinary  momentum  transfer  collisions  is  larger 
than  the  flow  scale  (3  meters). 

Optical  Data 

Two  optical  inctruments  were  used  in  this  experiment:  a  high  speed 
image  converter  camera  and  a  set  of  rotating  mirror  cameras.  Snapshots 
at  the  accelerator  taken  with  the  image  converter  camera  are  shown  in 
Fig.  4.  The  timing  and  exposure  of  the  camera  placed  at  the  rear  wall 
of  the  tank  facing  the  accelerator  were  adjusted  to  view  the  initial  stages  of 
breakdown.  The  exposure  time  of  50  nanoseconds  is  sufficiently  rapid  to 
resolve  asymmetries  on  the  scale  of  1  cm  (assuming  velocities  of  10?  cm/sec) 
From  the  photographs  in  Fig.  4,  it  appears  that  (1)  the  annmar  region  of 
the  accelerator  is  filled  with  luminous  plasma  as  early  as  \/Z\x  after  initia¬ 
tion  of  the  current  pulse  and  (2)  no  gross  (azimuthal)  asymmetries  are 
present  near  the  peak  of  the  current  cycle  (~  1  fj,  sec). 

In  addition,  measurements  were  made  of  the  convected  magnetic 
field  in  the  plasma  at  large  distances  from  the  source  (up  to  1  meter). 

These  measurements  were  made  at  various  azimuthal  locations  and  no 
gross  asymmetries  were  found  when  the  accelerator  was  loaded  with  four 
symmetrically  placed  wires. 
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Typical  data  from  the  ro'.-ting  mirror  cameras  are  shown  in  Fig.  5. 
These  are  x-t  diagrams  of  the  luminous  front  taken  at  a  background  pressure 
of  45  pt  with  and  without  an  applied  axial  field.  These  results  show,  with 
striking  clarity,  the  fundamental  result  of  this  experiment:  the  aluminum 
plasma  is  considerably  decelerated  when  a  magnetic  field  is  applied  albeit 
the  direction  of  the  field  is  parallel  to  its  motion.  It  must  be  emphasized 
that  this  phenomenon  is  observed  even  when  the  energy  of  the  plasma  is 
very  large  compared  to  the  energy  of  the  ambient  magnetic  field,  (see 
next  section).  This  means  that  the  ambient  magnetic  field  can  play  no 
direct  role  in  appreciably  decelerating  the  plasma,  but  can  only  act  as  an 
intermediary  in  the  momentam  transfer  process  resulting  from  charge 
exchange  and  direction  ionization  collisions. 

A  summary  of  the  mirror  camera  data  taken  at  45  /i  is  shown  in 
Fig.  6.  At  pressures  above  100/i,  the  aluminum  plasma  starts  to  pick  up 
the  ambient  gas  by  ordinary  momentum  transfer  collisions  and  therefore 
50ju  was  the  highest  pressure  for  which  the  effects  of  charge  exchange  and 
ionizing  collisions  alone  were  studied. 

The  data  in  Fig.  6  shows  position  versus  time  of  the  luminous  front 
for  four  different  applied  magnetic  field  intensities.  It  must  be  pointed  out 
that  the  aluminum  ion  plasma  docs  have  an  azimuthal  component  of  B  field 
which  appears  to  be  convected  with  it,  whether  or  not  there  is  an  applied 
axi'  1  bias  field.  Figure  7  shows  the  magnitude  of  this  field  as  a  function 
of  time  for  no  applied  field.  The  deceleration  of  the  aluminum  plasma 
(debris)  appears  to  be  a  function  of  the  magnitude  of  the  applied  axial  mag¬ 
netic  field  until  the  field  intensity  is  equal  to  200  gauss.  At  200  gauss  the 
gyro  radius  of  an  air  ion  is  either  50  cm  (at  (Pt  =14)  or  100  cm  (at  wt  =  28) 
when  the  ion  speed  is  equal  to  the  initial  debris  velocity.  The  experimental 
results  indicate  that  a  further  increase  in  B  field  or  flow  dimension  would 
not  affect  the  debris  motion. 

The  data  on  Fig.  7  for  Bz  =  0  shows  some  deceleration,  which  is 
possibly  due  to  two  effects. 

The  first  deals  with  the  manner  in  which  the  detector  (camera) 
records  the  position  of  the  leading  edge  of  the  luminous  front  and  this  is 
discussed  in  the  Appendix.  The  analysis  in  the  Appendix  is  for  an  ideal 
gas  expanding  in  a  vacuum.  It  is  shown  that  if  an  instrument  such  as  a 
camera  records  the  position  of  a  threshold  level  of  light  intensity  vs  time, 
its  trajectory  always  appears  to  decelerate. 

The  second  effect  causing  deceleration  is  that  even  with  no  applied 
field,  the  residual  magnetic  field  locked  in  the  plasma  primarily  in  the 
azimuthal  direction  and  therefore  perpendicular  to  the  motion.  Charge 
exchange  and  ionizing  collisions  can  couple  momentum  to  the  aluminum 
plasma  via  this  azimuthal  fiela  causing  the  plasma  to  decelerate  a  small 
amount  even  in  the  absence  of  an  applied  axial  field. 

The  deceleration  of  the  plasma  as  a  function  of  pressure  for  constant 
applied  axial  magnetic  field  (~  200  gauss)  is  shown  in  Fig.  8. 
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Fig.  8 


Aluminum  Plasma  Motion  at  Center  Line  vs  Time  with  Axial 
B  Field 


For  comparison,  x  vs  t  has  been  plotted  assuming  simple  momentum 
conservation  (snowplow  model  -  Eq.  (8)).  The  only  differences  between 
this  model  and  that  of  Eq.  (8)  is  that  the  neutral  gas  is  swept  up  by  charge 
exchange  and  ionizing  collisions  in  the  presence  of  a  magnetic  field  (instead 
of  ordinary  collisions).  To  obtain  the  best  fit  with  the  data,  a  cone  with 
"radius  to  height  ratio,  a  =  0.  6  was  used,  although  the  results  are  not  sensi¬ 
tive  to  the  choice.  Magnetic  field  data  presented  in  the  next  section  con¬ 
firms  this  geometry. 

Mirror  camera  data  was  also  obtained  at  15  /!  and  at  2  /Li .  These 
results  appear  to  scale  with  the  same  pressure  dependence  as  shown  in 

Fig.  8. 
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IV.  MAGNETIC  FIELD  DATA 


Magnetic  field  probes  were  used  first  to  confirm  the  velocities 
measured  by  the  mirror  cameras  and  second  to  determine  the  geometry 
and  magnitude  of  the  magnetic  field.  These  probes  were  placed  at  dif¬ 
ferent  axial  locations  to  measure  the  velocity  of  the  magnetic  disturbances. 
In  these  first  tests,  no  axial  field  was  applied,  so  the  only  field  present 
was  that  associated  with  the  plasma.  The  results  of  these  tests  can  be 
summarized  as  follows: 

1.  A  sharply  defined  {sec  Fig.  7),  pulse  was  observed  with 
no  radial,  Br  or  axial,  Bz  signals. 

2.  The  absolute  space-time  position  of  these  Bq  pulses  roughly 
correspond  to  the  mirror  camera  data. 

3.  The  time  difference  between  the  Bq  signal  for  various  axial 
stations  coincided  with  mirror  camera  data. 

-3/2 

4.  The  magnitude  of  the  Bq  field  scales  roughly  as  z  .  This 
fact  indicates  that  the  magnetic  energy  locked  in  the  aluminum  plasma  is 
constant  and  the  magnetic  energy  density  merely  decreases  as  the  volume 
of  the  bubble  increases. 

5.  The  magnetic  energy  associated  with  the  B$  field  is  many  orders 
of  magnitude  less  than  the  kinetic  energy  of  the  plasma.  The  value  for  B q 
is  usually  less  than  100  gauss,  which  when  evenly  distributed  over  a  volume 
30  cm  on  a  side  only  represents  a  total  magnetic  energy  of  (B^/2  p.  ) 
(volume)  -10”^  joules.  The  initial  kinetic  energy  of  the  aluminun?  plasma 
is  1/2  mvQ  —  2  x  10^  joules.  The  ratio  of  kinetic  to  magnetic  energy  = 

2  x  10^.  This  indicates  that  after  a  relatively  small  expansion,  the  ac¬ 
celerator  has  produced  a  plasma  with  only  kinetic  energy,  i.  e.  ,  there  is 
very  little  magnetic  field  pushing  the  plasma.  Consequently,  the  currents 
in  the  plasma  are  relatively  small  and  in  this  respect  the  plasm',  simulates 
nuclear  debris. 

The  remainder  of  the  magnetic  field  data  was  taken  with  an  applied 
axial  magnetic  field  (50-220  gauss)  in  order  to  measure  the  perturbation  in 
the  ambient  field.  Examples  of  the  B  field  data  taken  60  cm  from  the  ac¬ 
celerator  at  radii  of  10  cm  and  40  cm  are  shown  in  Figs,  9  and  10  and 
Fig.  11  respectively. 

Qualitatively,  the  dynamics  of  the  magnetic  field  bubble  are  as 
follows:  The  moving  plasma  pushes  the  ambient  axial  field  out  of  the  way 
and  creates  a  field  free  region.  After  the  plasma  passes,  the  field  closes 
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Fig.  10  Azimuthal  Magnetic  Field  Signals  (z  =  60  cm,  R  =  10  cn 
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behind.  A  composite  picture  of  the  magnetic  field  bubble  is  shewn  in 
Fig.  12.  The  bubble  is  approximately  40  cm  in  radius  and  this  agrees 
with  the  assumption  of  a  =*  0.  6  as  shown  in  Fig.  8.  This  picture  was 
drawn  from  data  taken  at  one  axial  location  (60  cm)  and  at  several  radial 
positions  and  "unfolded"  in  space  using  a  constant  velocity  (2.  5  cm/ji  see) 
obtained  from  mirror  camera  data. 

There  are  a  number  of  important  features  shown  in  Fig.  9  and  10. 
They  are: 

1.  The  rise  time  of  the  leading  edge  of  the  bubble  is  ~  8/i  sec 
which  at  a  velocity  of  2.  5  cm / \±  sec  corresponds  to  a  thickness  of  ~  20  cm. 
The  gyroradius  of  an  aluminum  ion  traveling  at  this  speed  in  a  field  of 
220  gauss  is  28  cm  and  thus  the  thickness  would  appear  to  be  comparable 
to  an  ion  gyro  radius.  If,  instead,  we  assume  that  at  60  cm,  the  aluminum 
plasma  has  a  large  amount  of  swept  up  air  within  it,  the  average  ion  mass 
could  be  lowered  (atomic  weight  of  dissociated  air  ~  14  as  compared  to 

~  27  for  aluminum). 

2.  The  ambient  axial  field  is  completely  excluded  for  approximately 

5  0  cm. 

3.  There  exists  a  spike  in  both  radial  and  axial  field  signals  which 
corresponds  exactly  to  the  time  that  the  azimuthal  field  which  is  recorded 
at  the  same  axial  location,  (see  Figs.  9  and  10)  has  a  sharp  increase.  The 
magnitude  of  the  azimuthal  field  is  approximately  equal  to  150  gauss  (when 
there  is  an  applied  axial  field)  while  the  spike  in  the  radial  and  axial  field 
traces  is  only  25  gauss. 

The  existence  of  the  disturbance  in  the  ambient  axial  field  for  a 
distance  approximately  equal  to  an  ion  gyro  radius  ahead  of  the  aluminum 
plasma  could  indicate  an  MHD  shock  (collisionless)  in  the  ambient  B 
field  and  plasma.  This  has  not  yet  been  verified  but  remains  as  an  in¬ 
teresting  point  for  further  study. 

4.  At  late  times  in  Fig.  9,  the  axial  field  trace  shows  slight 
enhancement  as  the  ambient  field  recloses  around  the  plasma. 

5.  In  Fig.  11,  at  a  larger  radial  position,  we  note  that  the  axial 
field  is  initially  enhanced  before  it  relaxes  back  to  its  original  state. 

This  effect  is  caused  by  the  bunching  up  of  the  axial  magnetic  field  lines 
at  the  edge  of  the  bubble  (see  Fig.  12).  We  note  that  in  this  high  velocity 
flow  (high  speed  w.  r.  t.  the  ambient  Alfven  speed)  the  wake  should  appear 
hypersonic  and  should  be  almost  flat  against  the  plasma  itself.  The  lines 
of  force  crowd  together  at  the  edge  and  the  radial  disturbance  does  not 
extend  far  from  the  plasma. 

Magnetic  field  data  at  this  same  pressure  (2  fi)  but  at  lower  values 
of  ambient  field  show  that  this  picture  is  unchanged.  At  higher  pressures, 
10/j.  -  45  ji,  it  is  found  that  the  plasma  still  excludes  the  ambient  field  over 
a  significant  volume  but  this  only  occurs  closer  to  the  accelerator,  i.e., 
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Fig.  12  Magnetic  Field  Profile 


before  the  plasma  is  greatly  decelerated.  At  the  observation  station  of 
z  =  60  cm,  R  =  10  cm,  75%  of  the  ambient  field  is  excluded  at  10  gi  and 
50%  at  45  n. 

There  is  also  a  qualitative  difference  in  the  high  pressure  data. 

The  data  is  less  reproducible  and  has  more  fluctuations  than  at  the  lower 
pressures.  It  is  possible  that  at  higher  pressures,  the  strong  rate  of  de¬ 
celeration  caused  instabilities  (e.  g.  ,  Taylor  instabilities)  which  might 
account  for  the  erratic  data.  It  is  interesting  to  point  out  that  this  experi¬ 
mental  facility  might  be  useful  in  the  simulation  of  late  time  debris  motion 
when  the  magnetic  field  is  collapsing  back  into  the  bubble. 


V. 


CONCLUSIONS 


This  experiment  has  demonstrated  that  a  magnetic  field  is  impor¬ 
tant  in  the  momentum  balance  of  the  moving  aluminum  plasma  even  though 
(1)  the  magnetic  field  energy  is  insignificant  when  compared  to  the  plasma 
kinetic  energy  and  (2)  the  ambient  magnetic  field  is  parallel  to  the  flow 
velocity.  The  way  in  which  the  magnetic  field  decelerates  the  plasma  is 
to  provide  a  momentum  coupling  for  charge  exchange  and  direct  ionization 
collisions. 

Several  aspects  of  high  altitude  nuclear  explosions  have  been  simu¬ 
lated  by  this  experiment  including  preionization  of  the  ambient  air,  a  high 
velocity  plasma  with  a  very  small  converted  magnetic  field  and  an  es¬ 
sentially  planar  front  parallel  to  the  ambient  field.  By  no  means  have  we 
attempted  to  simulate  the  entire  picture  of  high  altitude  nuclear  explosions, 
but  only  to  confirm  certain  assumptions  of  the  model  advanced  by  Workman.  ^ 
To  that  end,  we  have  confirmed  that  the  mo  ring  plasma  turns  the  ambient 
(axial)  field  until  it  lies  within  the  plane  of  the  bubble  and  that  a  field  free 
region  is  created.  Moreover,  the  streak  camera  photographs  show  that 
the  aluminum  plasma  gives  up  momentum  to  nearly  all  of  the  partially 
ionized  air  that  it  passes  over  when  the  ambient  axial  B  field  intensity 
just  large  enough  (200  gauss)  to  make  the  air  ion  gyro  radius  comparable 
to  the  flow  dimensions  (50-100  cm).  At  lower  B  field  intensities,  the  coupling 
appears  to  be  only  partial. 

There  are  many  details  of  the  flow  configuration  which  although  not 
essential  to  the  Workman's  model,  are  interesting  in  the  sense  of  plasma 
physic  research.  It  appears  from  the  magnetic  field  data,  that  the  ambient 
magnetic  field  is  turned  ahead  of  the  aluminum  plasma  region  and  that  the 
thickness  of  this  transition  region  is  of  the  order  of  an  ion  gyro  radius. 

Also,  the  magnetic  field  traces  show  structure  (e.  g,  ,  the  spike  of  azimuthal 
field  associated  with  the  plasma)  which  is  much  smaller  than  an  ion  gyro 
radius  indicating  that  perhaps  there  is  yet  a  smaller  relevant  scale  length. 
Finally  this  form  of  plasma  accelerator  has  proven  to  be  a  reliable,  well 
calibrated  tool  for  high  velocity  plasma  research  and  a  simple  method  has 
been  demonstrated  by  which  a  known  amount  of  material  can  be  vaporized, 
ionized  and  accelerated  reproducibly.  Since  only  one  half  the  capacitor 
bank  voltage  was  used  for  these  experiments,  it  should  be  possible  to  raise 
the  exhaust  velocity  to  perhaps  ~  2  x  10^  cm/sec  and  this  will  be  done  in 
the  near  future. 

Although  we  know  that  a  deceleration  mechanism  for  debris  along 
field  lines  does  exist,  the  second  prediction  of  Workman  concerning  the 
charge  exchange  leak  has  yet  to  be  investigated.  Confirmation  of  this 
charge  exchange  leak  is  one  of  our  remaining  goals. 
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LOW  PRESSURE  EXPANSION 


As  indicated  in  Fig.  6,  the  aluminum  plasma  decelerates  even  with 
no  applied  field.  In  Fig.  13,  data  is  shown  for  no  applied  field  at  15 /i  and 
45  fj,.  Aside  from  the  observation  that  both  curves  show  deceleration,  the 
lower  pressure  data  shows  a  higher  rate  of  deceleration  than  the  high  pres¬ 
sure  data.  These  two  phenomena  can  be  explained  in  a  qualitative  sense 
by  examining  the  characteristics  of  the  detector  (photographic  film). 

Film  requires  a  minimum  or  threshold  value  of  light  intensity,  Ij, 
to  visibly  record  an  event.  Ij  is  a  function  of  the  aluminum  plasma  density, 
n,  and  the  specific  functional  dependence  is  related  to  the  specific  radiation 
process  which  is  responsible  for  the  light  emission.  In  any  case,  1^  cor¬ 
responds  to  aluminum  density,  n^,  and  it  is  the  trajectory  of  n^  which  is 
shown  in  the  mirror  camera  photographs. 


As  a  model  for  the  calculation  of  this  trajectory,  consider  that  at 
time  t  =  0,  at  z  =  0,  the  aluminum  plasma  is  suddenly  given  the  one  di¬ 
mensional  velocity  distribution  as  shown  in  Fig.  14. 


f (v,  t=0,  z=0)  = 
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2  a 


1  /  \2 
2  (v"Vo} 


(Al) 


where  No  =  total  number  of  aluminum  particles 

a  =  v/kT  7 m  =  initial  sound  speed 
o  V  o  o 

v  =  z  component  of  velocity 

v  =  initial  velocity 
o 


Assuming  that  the  plasma  flows  into  a  field  free,  vacuum  region  and  that 
the  plasma  particles  do  not  interact  with  each  other,  the  density  distribution 
as  a  function  of  (z,  t)  will  be 
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Fig.  1 


Aluminum  Plasma  Motion  vs  Time  with  no  Applied  Field  and 
Pressure  as  a  Parameter 
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Fig.  14  Vacuum  Deceleration 
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A  three  dimensional  expansion  would  only  change  power  of  t  in  the  coef¬ 
ficient  of  the  exponential  from  one  to  three,  and  the  qualitative  results 
would  be  identical.  Denoting  by  z^,  the  point  which  corresponds  to  the 
threshold  density,  n^.  the  trajectory  of  this  point  is  given  by 


Zj(t)  =  v  t 
a  o 


1  +  m 
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2  log  T 
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where  m  =  v  /a 
o  o  o 
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v  o 


A  sketch  of  zj(t)  is  shown  in  Fig.  14  and  it  is  clear  from  the  curva¬ 
ture  that  the  acceleration  is  always  negative.  Furthermore,  if  the  vacuum 
region  is  replaced  by  a  low  pressure  background,  a  typical  pressure  de¬ 
pendence  for  n^  would  be 


n  ,  cc  pf- 
d 


p  <  0 


so  that 


log  T  / 1  cc  log  P 


It  can  then  be  shown  that  the  acceleration  of  z^t)  varies  as 
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dt 
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Thus,  on  the  basis  of  this  lov  pressure  expansion  into  a  given 
threshold  intensity  level  for  film,  (1)  zero  field  expansions  should  deceler¬ 
ate  more  than  into  those  of  higher  pressure. 

As  previously  stated,  there  is  always  a  residual  magnetic  field 
locked  into  the  plasma  so  that  this  vacuum  deceleration  mechanism  is 
never  singly  present  in  this  experiment.  However,  its  effect  in  this  experi 
ment,  as  well  as  in  any  optical  measurement  of  high  altitude  nuclear  ex¬ 
plosions,  must  be  recognized. 
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